 1:


 a := 1; evalf(a); b := 2; c := 3; E := a*z^2+b*z+c; niy := .4; G := E/(2*(1+niy)); Alfa11 := 1/E; Alfa12 := -(1*niy)/E; Alfa13 := -(1*niy)/E; Alfa22 := 1/E; Alfa23 := -(1*niy)/E; Alfa33 := 1/E; Alfa44 := 1/G; Alfa55 := 1/G; Alfa66 := 1/G; Liambda := 2*niy*G/(1-2*niy); Liambda11 := 2*G+Liambda; Liambda12 := Liambda; Liambda13 := Liambda; Liambda22 := 2*G+Liambda; Liambda23 := Liambda; Liambda33 := 2*G+Liambda; Liambda44 := G; Liambda55 := G; Liambda66 := G; A := Matrix(3, 3, [[Alfa11, Alfa12, Alfa13], [Alfa12, Alfa22, Alfa23], [Alfa13, Alfa23, Alfa33]]); B := evalm(1/A); Xalfa := -1/G; Yalfa := -1/G; Zalfa := -1/(2*G+Liambda); a1 := -(2*G+Liambda)/G; b1 := -1; c1 := -(G+Liambda)/G; d1 := -1-Liambda/G; a2 := -1; b2 := -(2*G+Liambda)/G; c2 := -(G+Liambda)/G; d2 := -1-Liambda/G; a3 := -(G+Liambda)/(2*G+Liambda); b3 := -(G+Liambda)/(2*G+Liambda); c3 := -G/(2*G+Liambda); d3 := G/(2*G+Liambda);

 2:


 beta13 := B[1, 3];
 beta23 := B[2, 3]; beta33 := B[3, 3];
 z := 0; evalf(beta13); evalf(beta23); evalf(beta33); evalf(Liambda44); evalf(Liambda55);
 z := .1; evalf(beta13); evalf(beta23); evalf(beta33); evalf(Liambda44); evalf(Liambda55);



